Critical Behaviour of One-particle Spectral Weights in the Transverse Ising Model 
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We investigate the critical behaviour of the spectral weight of a single quasiparticle, one of the key 
observables in experiment, for the particular case of the transverse Ising model. Series expansions 
are calculated for the linear chain and the square and simple cubic lattices. For the chain model, a 
conjectured exact result is discovered. For the square and simple cubic lattices, series analyses are 
used to estimate the critical exponents. The results agree with the general predictions of Sachdev. 
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There is currently much interest in quantum phase 
transitions in strongly correlated magnetic and/or elec- 
tronic systems in condensed matter physics. The physics 
in the vicinity of such quantum critical points can be de- 
scribed via general scaling theories, independent of any 
particular microscopic model. Sachdev [Jj has developed 
a scaling theory for the spectral weight (defined below), a 
quantity which can be measured experimentally, and thus 
provides a bridge between model calculations and real 
materials. Our goal in this letter is to verify Sachdev's 
predictions for a particular quantum spin model, the 
Ising model in a transverse field. To the best of our 
knowledge, this is the first such microscopic verification 
of the scaling form. 

A fundamental quantity to describe the dynamical 
properties of quantum magnets is the dynamical struc- 
ture factor 
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where the notation is standard, and the angular bracket 
denotes a thermal average or, at zero temperature, a 
ground state expectation value. As is well known, S(k, u>) 
can be related directly to the inelastic neutron scattering 
intensity, allowing a comparison between experiments on 
real materials and theoretical model calculations. 

Integrating over energy, we obtain the integrated 
or static structure factor, 
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which is just the Fourier transform of the spin-spin cor- 
relation function. 

At low temperatures the major contribution to S(k, ui) 
will often come from low energy quasiparticle excitations. 
It is then useful to write in the 'spectral form' (valid 
for a discrete spectrum) 
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where the sum is over all eigenstates with energies E\, 
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and the 'spectral weights' 5? (k) are given by 



sf (k) = nr(k)n£(k) 



with 
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We consider only the zero temperature case. 

Henceforth we consider only the 'one-particle' spectral 
weights 5i p (k), referring to a state of a single quasipar- 
ticle with momentum k. Near a quantum phase transi- 
tion [l| these quantities are expected to show universal 
scaling properties, which can be derived from a field the- 
oretic approach, independently of the details of a partic- 
ular model. The spectral weight can be written, up to 
constant factors, as 
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where A(k) is the 'quasiparticle residue', i.e. the residue 
of the Green's function at the quasiparticle pole, and 
E(k) is the quasiparticle energy at momentum k. Near 
the critical point, these functions are expected to scale 
as 



A(k) ~ ( 5c - 5 r, 

E(k c ) ~ (g c - 9 y 
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where g is a coupling or tuning parameter, k c is the crit- 
ical momentum at which the energy gap vanishes, and 
ry, v are the standard critical exponents. Then Si p is ex- 
pected to scale as 



Sip(k) ~ {go-gV, k^k c ; 
S lp (k c ) ~ (g e -g)to- 1)v . 
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Our aim is to confirm and demonstrate this scaling 
behaviour for the transverse field Ising model 
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where af = 2Sf are Pauli operators and the second sum 
is over nearest neighbour pairs. This model has a quan- 
tum phase transition at A = 1 in one dimension, and 
at some specific A c for each higher-dimensional lattice. 
Our approach is to derive series expansions in A for vari- 
ous one-particle spectral weights in the 'disordered' phase 
A < A c , via the linked-cluster techniques described by 
Oitmaa et al. and to analyse the series by standard 
methods to obtain the critical behaviour. 

The Linear Chain 
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FIG. 1: (Color online) The 1-particle spectral weight S^' p ~(k) 
for the linear chain model graphed against momentum k/n, 
for various couplings A. 
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FIG. 2: (Color online) As Fig. for the 1-particle spectral 
weight Sr p +(fe). 



FIG. 3: The integrated 1-particle structure factors $i and "3?2 
(corresponding to S±~, respectively) graphed as func- 

tions of coupling A. 



The transverse Ising chain model is exactly solvable, 
and expressions for the energy spectrum, magnetization, 
etc. have been given by Pfeuty Q|. We have computed 
series for the one-particle matrix elements H,^ p (k) for this 
model up to order A 22 . The leading terms for Q^ p (k) are 
given in Ref. 0; we can supply further coefficients on 
request. 

While analyzing these series for their critical be- 
haviour, we realized that they could be represented ex- 
actly, up to the order calculated, by the closed form ex- 
pressions: 

^) = ld-x^-^l (10) 

where A(fc) is just half the 1-particle energy E(k): 

A(k) = E(k)/2 = [1 + A 2 - 2Acos(fc)] 1/2 . (11) 

It is natural to conjecture that this result is exact to all 
orders. 

One can now write down any desired spectral weight 
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^ + p + ( fc ) = ^p'(fc) = ^(l-A 2 ) 1/4 (1 (14) 

These formulae give a very clear demonstration of the 
critical behaviour expected theoretically. The quasipar- 
ticle residues for S h , are respectively 

A(fc) = i(l-A 2 ) 1 /4(i T A(fc)) 2 (15) 

which vanish at the critical point A = 1 with exponent 
1/4 = f]v, as expected for the transverse Ising chain 
where rj — 1/4, v = 1. At the critical momentum k = 0, 
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the energy gap is 2(1 — A), and vanishes with exponent 
u=l. 

The behaviour of Si~(k) and S^~(k) for several A is 
shown in Figs. ^ and [21 We can see that both spec- 
tral weights slowly decrease towards zero as A — ► 1, ex- 
cept at the point k — where they both diverge, just 
as we should expect. Note that vanishes at 

A = 2 cos(fc). 

The integrated 1-particle structure factors are 
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where m = 4A/(1 + A) 2 , and K(m) and E{m) are the 
complete elliptic integrals of the first and second kinds, 
respectively 5]. Their behaviour is shown in Fig. |3 
where we can see that <J>2 has its maximum value at A = 
0.98512, and when A = 1, both $1 and $2 vanish, despite 
the fact that <Si p diverges at k = 0. 

From equations I|1U|) . one can also show 
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Note that whereas Sf*(k) diverges as {A — > 1, k 
Sfp(k) does not. 
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By Fourier transforming, we obtain the 1-particle con- 
tributions to the correlation functions (using equation 
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Now exact expressions for the 1-particle contributions 
to the correlation functions in the 2D classical Ising 
model have been obtained long ago by Wu, McCoy, Tracy 
and Barouch 0,0- Taking the anisotropic 'Hamiltonian 
limit' in which the classical Ising model reduces to 

the quantum Ising model ifTjl. one can show that the ex- 
pressions of Wu et al. reduce to equation (jH?)l (details 
will be given elsewhere). This confirms that the results 
we have given are in fact exact. 



(17) The Square Lattice 



The transverse Ising model in (2+1) dimensions lies in 
the universality class of the 3D classical Ising model, and 
so its critical exponents are expected to be rj = 0.0364(5) 
, v = 0.6301(4), from various estimates [Tfij . 

The square lattice is bipartite, and hence one finds 
that the Hamiltonian is symmetric under a spin rotation 
by 7r about the z axis on the B sublattice, followed by 
a coupling inversion A — > —A. Correspondingly, there 
are symmetrical critical points at couplings A = iA^, 
where A c has been estimated from series expansions X¥3 
at 0.32851(8), and from a finite-size scaling analysis [I^l 
at 0.32841(2). 

We have computed series for the 1-particle spectral 
weight <5fp (k) and the quasiparticle energy E{k) to order 
A 14 . The calculation involved a list of 4654 284 clusters, 
consisting of up to 15 sites. The leading terms are 



E(k) = 2 — 2A(cos k x + cos k y ) 
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Sip(k) = 1/4 + X(cosk x +cos k y )/4 + X 2 [12 cos(k x ) cos(k y ) + 3(cos(2 k x ) + cos(2 k y ))]/16 



-— [15(cos k x + cos k y ) + 30(cos(2fc :E ) cos k y + cosk x cos(2k y )) + 5(cos(3k x ) + cos(3fc a ))] + 0(A 4 ) (21) 



(further terms can be supplied on request). Note that 
because the lattice is bipartite, both quantities are sym- 
metric under the simultaneous transformations {k — > 
7r-k,A^-A}. 

The results of a standard Dlog Pade analysis for Sf£(k) 
at k = (0, 0) and k = (tt/2, tt/2) are shown in Tablefl At 



k = (0,0), where the energy gap vanishes, the table al- 
lows estimates A c = 0.3284(4) with exponent —0.605(5), 
compared to the expected exponent v(r\ — 1) = —0.607. 
At momentum k = (7r/2,7r/2), where the energy gap re- 
mains finite, we find A c = 0.34(3) with exponent 0.04(2) 
compared to the expected value vr\ = +0.0229. In the 
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TABLE I: Pole and residue of N/M Dlog Pade approximants to Sfp for the square lattice. Defective approximants are marked 
with an asterisk. 





N-2/N 


N-l/N 


N/N 


N+l/N 


N+2/N 








k = (0,0) 






N= 1 




4000^-0 8000 "1 


9703('-0 3659') 


3874Y-1 0759"; 


9847f-0 31 38"! 


N= 2 


0.2953(-0.4682) 


0.3263(-0.5900) 


0.3285(-0.6022) 


0.3281(-0.5988) 


0.3276(-0.5942) 


N= 3 


0.3286(-0.6028) 


0.3281(-0.5996) 


0.3319(-0.5789)* 


0.3278(-0.5969) 


0.3282(-0.6017) 


N= 4 


0.3238(-0.5113)* 


0.3279(-0.5978) 


0.3280(-0.5985) 


0.3288(-0.6122) 


0.3284(-0.6037) 


N= 5 


0.3280(-0.5985) 


0.3279(-0.5974)* 


0.3282(-0.6016) 


0.3283(-0.6035) 


0.3284(-0.6036) 


N= 6 


0.3283(-0.6019) 


0.3284(-0.6059) 


0.3284(-0.6053) 


0.3284(-0.6050) 




N= 7 


0.3284(-0.6053) 


0.3284(-0.6045) 














k = (7r/2,7r/2) 






N= 2 




0.6030(0.5455) 


0.6030(0.5455) 


0.3428(0.0570) 


0.3428(0.0570) 


N= 3 


0.6030(0.5455) 


0.6030(0.5455) 


0.3428(0.0570) 


0.3428(0.0570) 


0.3586(0.0747) 


N= 4 


0.4231(0.1781) 


0.3579(0.0735) 


0.3579(0.0735) 


0.3521(0.0654) 


0.3521(0.0654) 


N= 5 


0.3579(0.0735) 


0.3579(0.0735) 


0.3521(0.0654) 


0.3521(0.0654) 


0.3425(0.0495) 


N= 6 


0.3473(0.0576) 


0.3407(0.0459) 


0.3407(0.0459) 


0.3390(0.0424) 




N= 7 


0.3407(0.0459) 


0.3407(0.0459) 









latter case, the exponent does not agree particularly well 
with the expected value, but this inaccuracy might easily 
be attributed to the weakness of this singularity. Comple- 
mentary results follow for the antiferromagnetic singular- 
ity at A = —0.3284, with exponent 1^(77 — 1) at k = (7r,7r), 
and vq at k = (ir/2, ir/2). 

The Simple Cubic Lattice 

The transverse Ising model in (3+1) dimensions lies 
in the universality class of the 4D classical Ising model, 
where we expect the mean field exponents rj = , v = 
1/2, modulo logarithmic corrections 0. 

We have computed series for the 1-particle structure 
factor Sfp(k) and the quasiparticle energy E(k) to or- 
ders A 10 . The critical point has been obtained previously 
as A c = 0.19406(6). A Dlog Pade analysis of Sf£(k) 
at k = (0,0,0), where the energy gap vanishes, gives 



A c = 0.19406(8) with exponent -0.54(1), while for an- 
other momentum, where the energy gap remains finite, 
we find A c = 0.22(3) with exponent 0.03(2). Allowing for 
logarithmic corrections, these agree reasonably well with 
the expected values. 

In summary, then, we have shown that the one-particle 
spectral weights of the transverse Ising model satisfy the 
scaling behaviour at a quantum phase transition as pre- 
dicted by theory pj . For the linear chain we obtain new 
exact results; for the square and simple cubic lattices 
numerical estimates are obtained from series expansions. 
Further details will be published in a full length paper. 
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